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Abstract: Although the Suidae animals were one of the most prosperous mammals, their taxa and phylogeny were 
poorly studied. To bridge this research gap, the complete mitochondrial DNA cyt b sequences (1 140 bp)of one red river 
pig ( Potamochoerus porcus ), one bearded pig (Sus barbatus), and several Eurasian wild boars (Sus scorfa ) were deter- 
mined with direct PCR sequencing method. Incorporated with the sequences of other Suidae species obtained from Gen- 
Bank, the phylogentic trees of Suidae species were constructed by Neighbor-Joining and Most Parsimony method. Our re- 
sults showed that the babirousa was sister group of all other species. The African species clustered into a single clade, 
while the Eurasian species formed the other clade. The phylogenetic positions of the Sus species were consistent with their 
geographical distribution. Based on our phylogeographical results, the Sus species could be classified into three groups: 
(1) the relatively primitive animals from Philippine Islands (such as S. cebifrons and S. philippensis ), Sulawesi and it’s 
adjacent islands (S. celebensis), (2) the Eurasian wild boars (S. scrofa and the mysterious S. salvanius) , and (3) oth- 
er pigs (such as S. barbatus and S. verrucosus ), which might further diverged as the islands type (such as S. barbatus 
from Sumatra and Borneo) and the continental type (such as Malay Peninsula wild boar) . 
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The family Suidae animals were widely distributed in old world continents and Southeast Asian Islands. 
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Simpson (1945) recognized three subfamilies from the 
extant species: Suinae (including genus Sus , Potamo- 
choerus, and Hylochoerus ), Phacochoerinae ( Phaco- 
choerus ) and Babyrousinae ( Babyrousa ). The babirusa 
( Babyrousa babyrussa ) was distributed only on Sulawe- 
si and some neighboring islands in Indonesia, while 
Phacochoerinae (warthog) was distributed in tropical 
area of Africa. Among the Suinae, both Potamochoerus 
and Hylochoerus were only distributed in Africa. The 
Sus were one of the most widely distributed and pros- 
perous mammals, but the number of species in genus 
Sus is unclear. According to Groves (1981), at least 5 
species and 24 subspecies are recognized, while Hard- 
jasasmita (1987) recognized several additional species. 
Later, Groves (1997, 2001) and Oliver (2001) classi- 
fied the genus Sus into three species groups: (1) the 
S. verrucosus group, including the Javan warty pig S. 
verrucosus , the bearded pig S. barbatus (from Borneo, 
Sumatra, Bangka, the Riau Archipelago, Palawan and 
the Malay Peninsula), the Sulawesi pig S. celebensis , 


and the Visayan pig S. cebifrons; (2) the S. philip- 
pensis group, and (3) the S. scrofa group, including 
Eurasia wild boar S. scrofa and the pygmy hog S. sal- 
vanius . By analyzing morphologic and molecular data, 
Lucchini et al (2005) disapproved with this species 
group classification. In summary, the phylogenetic re- 
lationships among the species, species groups, genus 
and subfamiles of Suidae are not clear (Groves, 2001; 
Lucchini et al, 2005). To bridge this research gap, we 
analyzed the complete mtDNA cyt b sequences of the 
species in Suidae. 


1 Materials and Methods 
1.1 Sampling 


One red river pig Potamochoerus porcus, one 
Malaysia wild boar Sus barbatus , and several Eurasian 
wild boars and domestic pigs Sus scrofa , were collected 
for sequencing. The sequences of other species (Kijas 
& Andersson, 2001; Randi et al, 1996, 2002) were 
obtained from GenBank (Tab. 1). 


Tab. 1 Location and GenBank accession number of species used in the phylogenetic analysis 








Common name Species Location GenBank accession No. 
Peccaries Tayassu pecari America U66290 
Babirusa Babyrousa babyrussa Indonesia Z50106 
Babirusa B. babyrussa Indonesia AJ314559 
Babirusa B. babyrussa Indonesia AJ314560 
Babirusa B. babyrussa Indonesia AJ314545 
Babirusa B. babyrussa Indonesia AJ314546 
Warthog Phacochoerus africanus Africa: Tanzania AJ314547 
Warthog P. africanus Africa: Kenya AJ314548 
Warthog P. aethiopicus Africa: Kenya AJ314549 
Warthog P. aethiopicus Africa: Kenya AJ314550 
Warthog P. aethiopicus Africa: Kenya AJ314551 
Red river pig Potamochoerus porcus Africa DQ315602 
Javan warty pig Sus verrucosus AJ314552 
Javan warty pig S. verrucosus AJ314553 
Bearded pig S. barbatus AJ314554 
Bearded pig S. barbatus AJ314555 
Bearded pig S. barbatus AY534297 
Bearded pig S. barbatus Z50107 
Philippine warty Pig S. philippensis AY920905 
S. philippensis x S. barbatus AY920907 
Visayan warty pig S. cebifrons AY920906 
Bearded pig S. barbatus Malay Peninsula DQ315601 
Asian wild boar S. scrofa China DQ315597 
Asian wild boar S. scrofa China DQ315598 
Asian wild boar S. scrofa China DQ315599 
Asian domestic pigs S. scrofa China DQ315600 
Italian wild boar S. scrofa Italy AF304201 
Swedish wild boar S. scrofa Swedish AF304203 
Landrace S. scrofa domestic European NC-000845 
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1.2 DNA amplification and sequencing 

Genomic DNA was extracted by standard phenol/ 
chloroform method. The complete mtDNA cyt b se- 
quences were amplified and sequenced by the following 
primers L15280: (5'-GGA ATC TAA CCA CGA CCA 
A-3'), H15829: (5'-AAA AGC CCC CTC AGA TTC- 
3'); L15562: (5'-ATT ACG GAT GAG TTA TTC GCT 
ACC-3'), H16123 (5'-GGG TTT GCT GGG GTG TAG 
TTG TCT-3'); L15971 (5'-AAT CTC ATC AGA CAT 
AGA C-3'), H16543 (5'-AGT CTC CTT CCT TGA 
GTC TTG G-3'); L16254 (5'-GCC CAT ACT ACA 
CAC ATC CA-3'), H180 (5'-GTT AAT GCA CGA 
CGT ACA TAG-3'). L and H refer to light and heavy 
strands, respectively. PCR were performed in a 50 uL 
volume [500 mmol/L Tris-HCl (pH 8.3), 0.1% Triton 
X-100, 2.5 mmol/L KCl, 1.5 mmol/L MgCl, 0.5 
mmol/L of each dNTP, 10 pmmol/L of each primer, 
and 1 unit of Taq polymerase (Sagc) | following 35 cy- 
cles of Imin at 94 °C, 1 min at 58 °C, and 1 min at 72 
C. PCR products were purified on spin columns (Wat- 
son BioTechnologies Inc., Shanghai) and were directly 
sequenced for both strands by using BigDye™ Termina- 
tor Cycle Sequence Kit (ABI Applied Biosystems ) ac- 
cording to the manufacturer’s manual. 
1.3 Data analysis 

Sequences were edited and aligned by DNAstar 
package (DNASTAR Inc.). 


pecari, that the most closely related to Suidae, was 


The species, Tayassu 
used as the outgroup in the phylogenetic analysis. The 
neighbor-joining (NJ) tree (Saitou & Nei, 1987) was 
constructed using Kimura 2-parameter included in 
MEGA version 3.0 (Kumar et al, 2004). The most 
parsimony tree was constructed with PAUP 4.0 8 
(Swofford, 1998) using the branch-and-bound search 
with TBR branch-swapping option. Robustness of the 
nodes was assessed by the bootstrap method ( Felsen- 
stein, 1985) after 1 000 replications (bootstrap option 
with heuristic search in PAUP) by adding sequences 


random. 


2 Result and Analysis 


2.1 The phylogenetic relationships among the sub- 
families and genus 

To investigate the phylogenetic relationships among 
the subfamilies and genus, the phylogenetic trees of 
Suidae species were constructed by both Neighbor-Join- 
ing and Most Parsimony method. In NJ tree (Fig. la), 
Phacochoerus and Suidae clustered into a large branch, 
which became the sister group of the clade solely com- 
posed of Babyrousa. However, the sister relationship 
between Potamochoerus and Phacochoerus received low 


statistic support (BP = 28). MP analysis of the dataset 
(1 140 nucleotide characters with 363 variable sites and 
261 parsimony-informative sites) generated 10 MPTs 
(CI = 0.655, RI = 0.827, RC = 0.542, L = 667). 
The strict consensus tree of these 10 MPTs was pro- 
duced as shown in Fig. 2b. Both methods produced es- 
sentially identical tree topologies. In MP analysis, the 
sister relationship between Potamochoerus and Phaco- 
choerus received moderate statistic support (BP = 65). 
The phylogenetic status of Babyrousa was consistent 
with its special morphological characters, such as the 
more complex stomach (Langer, 1988) , the remarkable 
canines of the male (Groves, 1981), compared with 
other Suidae animals. Interestingly, a recently study 
(Bernhard & Lowden, 2004) on suid herpesviruses also 
showed that the virus Babyrousa babyrussa rhadinovirus 
1 (Bbab RHV-1), identified from Babyrousa, were 
more closely related to a ruminant gammaherpesvirus , 
bovine herpesvirus 4 (BoHV-4), than to the virus iden- 
tified from other suid animals, which were closely relat- 
ed to known beta- (porcine cytomegalovirus ) and 
gammaherpesviruses [porcine lymphotropic herpesvirus 
(PLHV) 1 and 3] of domestic pigs. So the molecular 
evidence, together with the morphological evidences 
support the suggestion that Babyroussa becomes an in- 
dependent family. In the large clade composed of Pha- 
cochoerus and Suinae, all Eurasian species were clus- 
tered into a single clade, while the African species 
formed another clade with moderate support. This 
makes the Sus paraphyletic, but seems in favor of the 
suggestions that the African species is monophyletic 
(Groves & Grubb, 1993) . 
2.2 The phylogenetic relationships among the Sus 
species 

The visayan wild pigs S. cebifrons branched off 
firstly from the Eurasian clade (Fig. 1), this was con- 
sistent with Lucchini’s (2005) results. Although the S. 
celebonse was not analyzed here, but it might have the 
similar phylogenetic status since it retains the same ple- 
siomorpic character with S. cebifrons (Groves, 1997). 
The next species diverged from the branch was the 
Philippine wild boar S. philippensis. The phylogenetic 
position of visayan and Philippine wild pigs showed that 
the species in Philippine Islands and sulawesi were 
more primitive than species from other places. The S. 
barbatus and S. verrucosus were clustered into a single 
clade, and became a sister group of the Eurasian wild 
boar S. scrofa which was clustered into another single 
clade . It’s noteworthy that this sister group relationships 
didn’t support the previous results that the S. barbatus 
and S. verrucosus were the ancestral population of 
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Fig. 1 Phylogenetic trees constructed by neighbor-joining method and most parsimony method 


The numbers on nodes are the bootstrap value (1 000 replicates). Both trees produce the identical topology. 


(a) NJ tree; (b) The strict consensus tree derived from the MP analysis. 


Eurasian wild boars (Larson et al, 2005). In the Sus 
barbatus clade, the pigs from Malay Peninsular occu- 
pied a distinct clade, while the animals from Islands 
clustered into another clade. 
2.3 The species groups classification and their 
phylogenetic relationships 

Groves (1997, 2001) classified the Sus species 
into three groups based on their morphological charac- 
ters. By analyzing both morphological and molecular 
data, Lucchini’s (2005) reclassified the species groups 
as follows: (1) the pigs from the Philippines (S. ceb- 
ifrons ) and Sulawesi (S. celebensis ), (2a) the western 
S. barbatus subspecies, (2b) a group consisting of S. 
scrofa , S. verrucosus , and S. cf. barbatus ahoenobar- 
bus . Our results showed that the phylogenetic relation- 
ships were consistent with their geographical distribu- 
tion, so the Sus species could be classified accordingly 
as the following three groups: (1) the relatively primi- 


tive animals from Philippine Islands (such as S. ceb- 


ifrons and S. philippensis ), Sulawesi and it’s adjacent 
islands ( S. celebensis), (2) the Eurasian wild boars 
(S. scrofa and the mysterious S. salvanius ), and (3) 
other pigs (such as S. barbatus and S. verrucosus ) , 
which might further diverged as the islands type (such 
as S. barbatus from Sumatra and Borneo) and the con- 
tinental type (such as Malay Peninsula wild boar). 
Each species group was similar but different from previ- 
ous classification (Groves, 1997, 2001; Lucchini’s, 
2005) . 

In general, our phylogenetic analysis showed that 
the Africa species, the Sulawesi Babyrousa, and the 
Eurasia species, each formed an independent group 
that was consistent with their geographical distribution . 
Furthmore, the phylogenetic relationships among the 
species in Suidae were also consistent with their geo- 
graphical distribution. This finding would be beneficial 
for further research on their early dispersal history, tax- 


onomy and conservation. 
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